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ABSTRACT 

We present 14 N° measurements from our HIRES/Keck database of damped Lya abundances. These 
data are combined with measurements from the recent and past literature to build an homogeneous, 
uniform set of observations. We examine photoionization diagnostics like Fe ++ and Ar° in the majority 
0^ ■ of the complete sample and assess the impact of ionization corrections onN/a and a/H values derived 

from observed ionic column densities of N°, Si + , H°, and S + . Our final sample of 19 N/a, a/H pairs 
appears bimodal; the majority of systems show N/a values consistent with metal-poor emission regions 
in the local universe but a small sub-sample exhibit significantly lower N/a ratios. Contrary to previous 
studies of N/a in the damped systems, our sample shows little scatter within each sub-sample. We 
consider various scenarios to explain the presence of the low N/a sightlines and account for the apparent 
i-C ' bimodality. We favor a model where at least some galaxies undergo an initial burst of star formation 

with suppressed formation of intermediate-mass stars. We found a power-law IMF with slope 0.10 or 
a mass cut of ~ 5 — 8M Q would successfully reproduce the observed LN-DLA values. If the bimodal 
distribution is confirmed by a larger sample of measurements, this may present the first observational 
c/3 , evidence for a top heavy initial mass function in some early stellar populations. 



Subject headings: galaxies: abundances — galaxies: chemical evolution — quasars : absorption lines 
nucleosynthesis 
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A complementary approach toward tracing the detailed 
enrichment history of protogalaxies is to consider the rel- 
ative abundance of nitrogen. Nitrogen is mainly produced 
in the six steps of the CN branch of the CNO cycles within 
H burning stellar zones, where 12 C serves as the reac- 
tion catalyst (see a textbook like Clayton 1983 or Cowley 
1995 for a nucleosynthesis review). Three reactions occur 
to transform 12 C to 14 N: 12 C(p,7) 13 N(/3+^) 13 C(p, 7 ) 14 N, 
while the next step, 14 N(p,7) 15 0, depletes nitrogen and 
has a relatively low cross-section. The final two reactions 
in the cycle transform 15 O to 12 C. Since the fourth reac- 
tion runs much slower than the others, the cycle achieves 
equilibrium only when 14 N accumulates to high levels, and 
so one effect of the CN cycle is to convert 12 C to 14 N. 

One issue in nitrogen evolution is to discover the source 
of the carbon which is converted into nitrogen, and of 
any oxygen which can contribute through the (slow) side 
chain 16 0(p,7) 17 F(/3+,^) 17 0(p,a) 14 N. For example, stars 
may produce their own carbon (and some oxygen) during 
helium burning, and the carbon (and perhaps oxygen) is 
subsequently processed into 14 N via the CN(O) cycle. In 
this case, nitrogen production is independent of the initial 
composition of the star in which it is synthesized and is 
referred to as primary nitrogen. On the other hand, stars 
beyond the first generation in a galactic system already 
contain some carbon and oxygen, inherited from the in- 
terstellar medium out of which they formed. The amount 
of nitrogen formed from CNO cycling of this material will 
then be proportional to its C abundance (and also its O 
abundance, if the CNO cycling proceeds long enough to 
deplete the oxygen) and is known as secondary nitrogen. 
To zeroth order, then, primary nitrogen production is in- 
dependent of metallicity, while secondary production is lin- 
early proportional to metallicity. The effects of metallicity 
on nitrogen production is one facet of nitrogen evolution 
we would like to understand clearly. 

A second issue of concern in the origin of nitrogen is to 
identify the portion of the stellar mass spectrum which is 
most responsible for 14 N production. While only massive 
stars (M>8M Q ) can produce the necessary internal tem- 
peratures required to synthesize most heavy elements, in 
the case of nitrogen intermediate-mass stars (IMS; 1-8M Q ) 
are also massive enough to attain required temperatures 
for its synthesis. In fact, both observational evidence and 
theoretical predictions strongly indicate that significant 
amounts of nitrogen are produced in intermediate-mass 
stars (M<8M Q ), while some is also synthesized in massive 
stars (Maeder 1992; Woosley & Weaver 1995; van den Hock 
& Groenewegen 1997; Henry, Kwitter, & Bates 2000; Henry, 
Edmunds, & Koppen 2000, hereafter HEK00; Marigo 2001; 
Meynet & Maeder 2002; Siess, Livio, & Lattanzio 2002). 
And since the timescale for IMS evolution is longer than 
for massive stars, nitrogen production and ejection into 
the interstellar medium may be significantly delayed with 
respect to oxygen and other heavy elements that are pro- 
duced primarily in massive stars. This in turn means that 
nitrogen evolution is affected by the form of the initial 
mass function (IMF), time dependence of the star forma- 
tion rate, and the effective stellar nitrogen yields of both 
massive and intermediate-mass stars. 

Local measurements of emission line regions have exam- 
ined the abundance of nitrogen relative to O and other a- 



elements as a function of metallicity (see Henry & Worthey 
2000 for a compilation). These observations have provided 
evidence for the contributions of both the primary and sec- 
ondary nitrogen mechanisms described above to the cosmic 
abundance evolution of this clement. For example, one ob- 
serves (i) a 'plateau' of N/a measurements at [a/H] < — 1 
consistent with primary nitrogen formation, i.e., indepen- 
dent of metallicity; and (ii) a rise in N/a with increasing 
metallicity above [a/H] « — 1 suggestive of secondary ni- 
trogen formation. Similar observations can be performed 
with the damped Lya systems (Pettini et al. 1995; Lu et 
al. 1998, hereafter L98; Centurion et al. 1998). Because 
damped systems tend to have low metallicity (Z < Zq/10) 
their N/a measurements probe the primary regime of ni- 
trogen production. In addition, since these systems are 
relatively young and may not have achieved a stage of 
steady evolution because of differences in stellar evolution 
time scales over the stellar mass spectrum, N/a may al- 
low us to gauge the effects of time delay and the role of 
IMS in nitrogen production. And unlike the a/Fe ratios 
described above, the N/a observations are largely free of 
the uncertainties due to differential depletion because of 
the mild refractory nature of N, O, S, and Si. In fact the 
most significant source of uncertainty is due to the effects 
of photoionization: measurements of ions N°, O , S + , H°, 
and Si + must be converted to elemental abundances. 

The initial studies of N/a in damped systems suggested 
a significant scatter in the N/a ratios (L98) with a few 
systems showing values below the plateau defined by low 
metallicity H II regions. These observations are difficult 
to interpret as either purely primary or secondary nitrogen 
formation and a combination of the two channels was sug- 
gested (L98). In this paper, we will reexamine these con- 
clusions. We present 14 nitrogen measurements from our 
UCSD HIRES/Keck I damped Lya database (Prochaska 
et al. 2001; hereafter, P01). and combine these data with 
six from previous studies and an additional five measure- 
ments from the recent literature. We restrict our discus- 
sion to high resolution echelle observations because (1) 
the NI transitions lie within the Lya forest where line- 
blending is important; (2) several NI transitions are sub- 
ject to line saturation; and (3) the NI 1134 and 1200 
triplets are very closely spaced. We limit the analysis to 
the damped Lya systems where the effects of photoioniza- 
tion are more likely to be small. We synthesize the abun- 
dance measurements from these various studies, investi- 
gate the effects of photoionization, and present a homoge- 
neous, uniform sample of N/a and a/H measurements. We 
then consider various models of N production and describe 
the implications for the star formation histories within the 
damped Lya systems. 

2. THE UCSD HIRES DATABASE 

All of the systems comprising the UCSD database were 
observed on the Keck I telescope with the HIRES spec- 
trograph (Vogt et al. 1994). Most of the observations are 
described in P01 and we refer the reader to that paper 
for more details on the acquisition and processing of the 
data. For several of the systems presented in this section, 
we also rely upon measurements made with the UVES 
spectrograph (Dekker et al. 2000) on the VLT. In the fol- 
lowing subsections, we present the nitrogen transitions as 
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well as other transitions relevant to our analysis. These 
include intermediate-ion transitions like FeUI 1122 which 
help reveal the ionization state of the gas and a-element 
transitions like Sill 1808 which are essential for assess- 
ing the enrichment history. The following sub-sections 
present the figures and ionic column density measurements 
for each system. For ions with multiple transitions mea- 
sured, we report the variance-weighted mean. Table 1 lists 
the atomic data adopted throughout the paper where col- 
umn 2 lists the rest wavelength, column 3 gives the oscil- 
lator strength, and column 4 lists the reference. Complete 
abundance measurements made from the new observations 
will be provided in our abundance database 4 . 



Table 1 
ATOMIC DATA 



Transition 


A 


f 


Rcf 


OI 976 


Q7fi 4481 


00330000 


I 


NTTT 989 


989.7990 


10660000 

KJ . ± KJ UUUuU \) 


I 


Sill 989 

Ulll cue 


989.8731 


1 3300000 


1 


STTT 1fl12 


1 01 2 5020 


03550000 

KJ • KJiJ\J\jKJ\JKJKJ 


I 


OI 1039 


1039.2304 


0.00919700 


1 


Arl 1048 


1048.2199 


0.26280000 


3 


Fell 1063 


1063.1760 


0.06000000 


1 


Fell 1064 


1063.9718 


0.00371800 


2 


Arl 1066 


1066.6600 


0.06747000 


3 


Nil 1083 


1083.9900 


0.10310000 


1 


Fell 1121 


1121.9748 


0.02020000 


2 


FeUI 1122 


1122.5260 


0.16200000 


3 


NI 1134 


1134.1653 


0.01342000 


1 


NI 1134 


1134.4149 


0.02683000 


1 


NI 1134 


1134.9803 


0.04023000 


1 


Fell 1143 


1143.2260 


0.01770000 


2 


Fell 1144 


1144.9379 


0.10600000 


2 


SIII 1190 


1190.2080 


0.02217000 


1 


Sill 1190 


1190.4158 


0.25020000 


1 


Sill 1193 


1193.2897 


0.49910000 


1 


NI 1199 


1199.5496 


0.13000000 


1 


NI 1200 


1200.2233 


0.08620000 


1 


NI 1200 


1200.7098 


0.04300000 


1 


Silll 1206 


1206.5000 


1.66000000 


1 


HI-A 1215 


1215.6701 


0.41640000 


1 


SII 1250 


1250.5840 


0.00545300 


1 


SII 1253 


1253.8110 


0.01088000 


1 


SII 1259 


1259.5190 


0.01624000 


1 


OI 1302 


1302.1685 


0.04887000 


1 


Sill 1304 


1304.3702 


0.09400000 


4 


ALU 1670 


1670.7874 


1.88000000 


1 


Sill 1808 


1808.0130 


0.00218600 


11 


A1III 1854 


1854.7164 


0.53900000 


1 


A1III 1862 


1862.7895 


0.26800000 


1 



References. — Key to References - 1: Morton (1991); 
2: Howk et al. (2000); 3: Morton (2002); 
4: Tripp ct al. (1996); 11: Bcrgeson & Lawlcr (1993) 



With few exceptions, we calculate ionic column densities 
with the same procedures as described in P01, i.e., we im- 
plement the apparent optical depth method (Savage and 
Sembach 1991) for the majority of measurements. This 
technique, however, cannot be applied to transitions where 
line-blending or saturation is significant. This occurs more 
frequently for transitions presented in this paper because 
the majority arise within the Lya forest. For these cases, 
we either present upper limits on the measured column 

4 http://kingpin.ucsd.edu/~hiresdla 



densities or have performed a full line-profile analysis with 
the VPFIT software package kindly provided by R. Car- 
swell and J. Webb. Note that the tables report [X/H] 
measurements for the low-ion species having adopted the 
meteoritic solar abundances from Grevesse et al. (1996) 
except for N and Ar where we adopt their photometric val- 
ues and O where we adopt the revised value from Holweger 
(2001). 

2.1. PH957, z = 2.309 

Although this system is included in the UCSD/HIRES 
database, our observations do not cover the NI transi- 
tions. Therefore, we adopt the measurements of N° and 
several other ions from the UVES observations studied 
in Dessauges-Zavadsky, Prochaska, & D'Odorico (2002b): 
logA^N ) = 15.03±0.02,logAT(HI) = 21.37±0.08, log AT(S+) = 
15.11 ±0.01. 



Table 2 

IONIC COLUMN DENSITIES: Q0201+36, z = 2.463 



Ion 


A 


AODM 


Nadopt 


[X/H] 


HI 


1215.7 


20.380 ± 0.045 






N I 


1134.1 




> 15.000 


> -1.350 


N II 


1083.9 


> 15.303 






Si II 






15.532 ±0.010 


-0.408 ± 0.046 


S II 


1253.8 


15.287 ±0.011 


15.287 ±0.011 


-0.293 ±0.046 


Ar I 


1048.2 


14.077 ± 0.030 


14.077 ± 0.030 


-0.823 ± 0.054 


Fe II 






15.008 ± 0.004 


-0.872 ±0.045 


Fc III 


1122.5 


13.855 ±0.023 







2.2. Q0201+36, z = 2463 

We present a series of metal-line profiles for this damped 
system in Figure 1 including NI, FeUI, and Arl. Because 
the velocity profiles extend over 200 km/s, the two NI 
triplets are significantly self-blended. To estimate the N° 
column density, we have performed a line-profile analy- 
sis of these transitions with the VPFIT software pack- 
age. Unfortunately, the blending is so severe and com- 
plicated that we could not achieve a satisfactory solution. 
In turn, we can only place a very conservative lower limit 
of AT(N°) > 10 15 cur 2 . Note, the N(Fe++) value listed in 
Table 2 refers to the velocity region —120 km s _1 < v < 
40 km s _1 which does not include the strong feature at 
v « —150 km s" 1 , a probable line-blend. Similarly, we re- 
port an N(Fe + ) value corresponding to the same velocity 
region to enable an evaluation of the Fe ++ /Fe+ ratio (§ 3). 

2.3. Q0336-01, z =3.062 

Figure 2 presents the NI transitions observed for the 
damped Lya system at z = 3.062 toward Q0336-01. Un- 
fortunately all of the N I profiles are at least partially sat- 
urated and we can only set a lower limit on the N° column 
density. Table 3 lists the ionic column density measure- 
ments and the variance- weighted means. From FcIII 1122, 
we report an upper limit to N (Fe ++ ) because of the pres- 
ence of a significant line blend at v > km s _1 . Similarly, 
we classify our N(N + ) measurement as an upper limit be- 
cause we suspect line-blending with the Lya forest. 
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Fig. 1. — Velocity plot of NI and intermediate-ion transitions for the damped Lyo system at z = 2.463 toward Q0201+36. For comparison, 
we also plot the Si II 1808 profile. The vertical line at v = corresponds to z = 2.4628. 
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Fig. 2. — Velocity plot of the metal-line transitions for the damped Lya system at z = 3.062 toward Q0336-01. The vertical line at v = 
corresponds to z = 3.062078. 



6 



THE UCSD HIRES/KECK I DAMPED Lja ABUNDANCE DATABASE IV. 




-100 100 -100 

Relative Velocity (km s -1 ) 



Fig. 3. — Velocity plot of the metal-line transitions for the damped Lya system at z = 3.025 toward Q0347-38. The vertical line at v = 
corresponds to z = 3.0247. 
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Table 3 

IONIC COLUMN DENSITIES: Q0336-01, z = 3.062 



ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


21.200 ±0.100 






N I 


1134.9 


> 15.037 


> 15.037 


> -2.133 


N I 


1200.7 


> 15.014 






N II 


1083.9 


< 14.371 






S II 


1250.5 


15.118 ±0.022 


14.994 ± 0.011 


-1.406 ± 0.101 


S II 


1259.5 


14.970 ± 0.012 






S III 


1012.5 


< 14.650 


14.994 ± 0.011 


-1.406 ± 0.101 


Ar I 


1048.2 


> 13.939 


> 13.939 


> -1.781 


Fe II 






14.905 ± 0.033 


-1.795 ± 0.105 


Fc III 


1122.5 


< 14.238 










Table 4 




IONIC COLUMN DENSITIES: Q0347-38, 


z = 3.025 


Ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


20.626 ± 0.005 






C III 


977.0 


> 14.385 






N I 






14.890 ± 0.031 


-1.706 ±0.031 


Si II 


1808.0 


15.016 ±0.026 


15.016 ± 0.026 


-1.170 ± 0.026 


Si III 


1206.5 


> 13.447 






S II 






< 14.760 


< -1.066 


S III 


1012.5 


13.849 ±0.123 






Ar I 


1048.2 


> 14.031 


14.276 ± 0.035 


-0.870 ± 0.035 


Ar I 


1066.6 


14.276 ± 0.035 






Fe II 






14.503 ± 0.007 


-1.623 ±0.009 


Fc III 


1122.5 


13.140 ±0.080 







2.4. Q0347-38, z = 3.025 

The ionic transitions for the damped system toward 
Q0347-38 are presented in Figure 3. This damped sys- 
tem has been carefully analysed by several authors in- 
cluding Levshakov et al. (2002) who analysed the high 
quality UVES/VLT commissioning data of this quasar. 
In the following, we will adopt their HI measurement of 
log iV(IiI) = 20.626 ±0.005 and implement our own metal- 
line analysis except for nitrogen. The UVES observa- 
tions include many more NI transitions than our own 
and have significantly higher S/N ratio. We stress that 
the 7V(Si+) value that we have measured is « 0.25 dex 
lower than the value listed in Levshakov ct al. (2002). 
We believe the difference is due to a significant telluric 
line-blend in the UVES spectrum. Finally, the Fe ++ col- 
umn density reported in Table 4 refers to the velocity re- 
gion —20 km s _1 < v < 20 km s _1 and may be considered 
an upper limit because of the likelihood of blending with 
neighboring, coincident absorption features. 
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Fig. 4. — Velocity plot of the metal-line transitions for the damped 
Lya system at z = 3.017 toward HS0741+47. The vertical line at 
v = corresponds to z = 3.017399. 



Table 5 

IONIC COLUMN DENSITIES: HS0741+4741, z = 3.017 



Ion 


A 


AODM 


-'Vadopt 


[X/Hj 


HI 


1215.7 


20.480 ±0.100 






N I 


1134.4 


14.158 ±0.015 


13.977 ±0.010 


-2.473 ±0.100 


N I 


1134.9 


13.926 ±0.017 






N I 


1200.7 


14.012 ±0.012 






N II 


1083.9 


< 14.498 






Si II 


1020.6 


14.162 ±0.019 


14.354 ± 0.003 


-1.686 ±0.100 


Ar I 


1048.2 


13.134 ±0.020 


13.134 ±0.020 


-1.866 ±0.102 


Fe II 






14.052 ± 0.005 


-1.928 ±0.100 


Fc III 


1122.5 


< 13.290 







2.5. HS0741+47, z = 3.017 

The metal-linc profiles for this damped system are shown 
in Figure 4 which includes NI, Sill, Felll, Nil, and Arl 
transitions. Because we expect the Fe III 1122 and N II 1083 
profiles are dominated by blends with coincident Lya clouds, 
wc report their column densities as upper limits (Table 5). 
These limits were derived by integrating the velocity pro- 
files over the velocity intervals —20 km s" 1 < v < 30 km s" 1 
and —40 km s" 1 < v < 20 km s" 1 respectively. 

2.6. Q0930+28, z = 3.235 



The nitrogen abundance of this damped system was pre- 
viously analysed by L98. In their paper, they report an HI 
column density of 10 20 18 cm" 2 based on their observations 
of the Lya transition. We find the Lya profile is best fit 
by two HI absorbers at two redshifts: (1) z = 3.234946, 
iV(HI) = 10 20 ' 3 , and (2) z = 3.246338, N(HI) = 10 20 2 . 
We estimate a 0.1 dex error for each of these HI mea- 
surements. In the following, we restrict the analysis to 
the z = 3.234946 damped system because we have signif- 
icant concerns on the ionization state of the 'sub-DLA'. 
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Fig. 5. — Velocity plot of the metal-line transitions for the damped 
Lya system at z = 3.235 toward Q0930+28. The vertical line at 
v = corresponds to z = 3.2353. 
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Fig. 6. — Velocity plot of the metal-line transitions for the damped 
Lya system at z = 2.466 toward Q1223+17. The vertical line at 
v = corresponds to z = 2.466083. 




Table 6 

IONIC COLUMN DENSITIES: Q0930+28, z = 3.235 



Ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


20.300 ±0.100 






N I 


1134.1 


13.951 ±0.057 


13.740 ± 0.013 


-2.530 ±0.101 


N I 


1199.5 


13.659 ±0.022 






N I 


1200.2 


13.762 ±0.022 






N I 


1200.7 


13.884 ± 0.024 






N II 


1083.9 


13.655 ± 0.033 






Si II 






13.888 ± 0.021 


-1.972 ± 0.102 


Si III 


1206.5 


> 13.481 






S III 


1012.5 


13.837 ±0.089 






Ar I 


1066.6 


< 12.960 


< 12.960 


< -1.860 


Fe II 


1121.9 


13.549 ± 0.081 


13.695 ± 0.015 


-2.105 ±0.101 


Fe II 


1144.9 


13.767 ±0.032 






Fe III 


1122.5 


13.141 ±0.032 







Figure 5 shows the intermediate- ion and Arl transitions 
for this damped system and Table 6 lists the ionic column 
densities from our and the L98 analysis. 

2.7. Q1223+17, z = 2466 

We reveal the NI, Sill, and Felll transitions for the 
damped system toward Q1223+17 in Figure 6. The pro- 
files for this system are spread over 100 km/s such that the 
members of the NI 1134 triplet are blended with one an- 
other. We first estimated the NI column density from this 



Table 7 

IONIC COLUMN DENSITIES: Q1223±17, z = 2.466 



Ion 



A 



AODM 



HI 
N I 
N II 
Si II 
Fe II 

Fe III 1122.5 



1215.7 
1083.9 



21.500 ±0.100 
> 14.794 

< 14.082 



~WJW 



14.830 ±0.179 

15.467 ±0.008 
15.157 ±0.022 



-2.640 ±0.205 

-1.593 ±0.100 
-1.843 ±0.102 



triplet by integrating the apparent optical depth over all 
3 transitions and adopting a summed oscillator strength. 
We then performed a detailed line-profile analysis of the 
lines with the VPFIT software package. Unfortunately, at 
a S/N ratio of 5 we could not constrain the N(N°) value 
to better than 0.2 dex. For this reason, we have chosen to 
remove it from our analysis. 



2.8. 



-17, 



1.776 



Q1331- 

Although we adopt the Si + column density presented 
in P01, we rely on the analysis of Dessauges-Zavadsky, 
Prochaska, & D'Odorico (2002b) for the majority of mea- 
surements considered in this paper. Their UVES-VLT ob- 
servations cover the NI transitions and intermediate- ions 
which lie blueward of our wavelength coverage. 
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Fig. 7. — Velocity plot of the metal-line transitions for the damped Ly« system at z = 2.827 toward Q1425+60. The vertical line at v = 
corresponds to z = 2.8268. 
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Fig. 8. — Velocity plot of the metal-line transitions for the damped Ly« system at z = 2.844 toward Q1946+76. The vertical line at v = 
corresponds to z = 2.8443. 
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Table 8 

IONIC COLUMN DENSITIES: Q1425+6039, z = 2.827 



ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


20.300 ± 0.040 






N I 


1134.4 


14.707 ±0.014 


14.707 ±0.014 


-1.563 ± 0.042 


N I 


1200.2 


> 14.623 






N II 


1083.9 


> 14.715 






Si II 


1190.4 


> 14.751 


> 14.826 


> -1.034 


Si II 


1526.7 


> 14.826 






S III 


1012.5 


14.157 ±0.108 






Ar I 


1066.6 


< 13.426 


< 13.426 


< -1.394 


Fe II 


1143.2 


14.255 ± 0.039 


14.471 ± 0.006 


-1.329 ± 0.040 


Fc III 


1122.5 


< 14.018 







2.9. Q1425+60, z = 2.827 

L98 presented a lower limit on the N° column den- 
sity based on their HIRES spectra of the NI 1200 triplet. 
Our independent HIRES observations (Figure 7) provide 
a measurement of N° which includes the NI 1134 triplet 
as listed in Table 8. Our adopted value is just below their 
limit and we expect the difference is due to blending in 
the Lya forest. We note that the stronger N I transitions 
exhibit a small but significant feature at v w —90 km s _1 
which is evident in the Sill profiles as well as the Nil 1083 
profile. On the other hand, we suspect that the Felll 1122 
profile is significantly blended with a coincident Lya cloud 
and place an upper limit to iV(Fe ++ ). Unfortunately, both 
the L98 and our observations only provide a lower limit to 
Si + which serves as the a-element for this system. 

2.10. GB1759+ 75, z = 2.625 

A detailed analysis of this system, including NI, was 
presented in (Prochaska et al. 2002, hereafter P02). We 
refer the reader to that paper for figures and tables. 



Table 9 

IONIC COLUMN DENSITIES: Q1946±7658, z = 2.844 



Ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


20.270 ± 0.060 






N I 


1134.1 


< 13.214 


12.588 ± 0.039 


-3.652 ± 0.072 


N I 


1134.4 


< 12.810 






N I 


1134.9 


12.287 ±0.141 






N I 


1199.5 


< 12.860 






N I 


1200.2 


12.673 ± 0.038 






N II 


1083.9 


13.568 ±0.005 






N III 


989.7 


12.688 ± 0.063 






O I 






14.820 ± 0.007 


-2.320 ±0.060 


O I 


1039.2 


14.813 ±0.008 


14.820 ± 0.007 


-2.320 ± 0.060 


Si II 






13.604 ± 0.005 


-2.226 ± 0.060 


Si III 


1206.5 


12.860 ± 0.005 






S III 


1012.5 


< 13.476 






Fe II 






13.242 ± 0.009 


-2.528 ± 0.061 


Fe III 


1122.5 


< 13.216 







2.11. Q1946+ 76, z = 2.844 

The nitrogen abundance of this system has been ex- 
amined in several papers (Lu et al. 1996; L98). Our in- 
dependent observations (Kirkman & Tytler 1997) imply 
a value in good agreement with the best of the previous 



measurements: log7V(N°) = 12.58 ±0.04. In Figure 8, we 
present the N I transitions as well as several other metal- 
lines relevant to the system's ionization state and Table 9 
summarizes the column density measurements. We re- 
port a value for N(N + ) based on integrating the blended 
Nil 1083 profile over the velocity region — lOkms -1 < 
v < 10 kms -1 . Finally, the reported Si ++ value was de- 
rived by integrating the Si III 1206 profile over the velocity 
interval —10 km s _1 < v < 10 km s _1 and might be con- 
sidered an upper limit given the likelihood of line-blending. 
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Fig. 9. — Velocity plot of the metal-line transitions for the damped 
Lya system at z = 2.538 toward Q2344+12. The vertical line at 
v = corresponds to z = 2.53790. 



Table 10 

IONIC COLUMN DENSITIES: Q2344+12, z 



2.538 



Ion 


A 


AODM 


Nadopt 


[X/Hj 


HI 


1215.7 


20.360 ±0.100 






N I 


1134.9 


13.742 ±0.088 


13.779 ±0.029 


-2.551 ±0.104 


N I 


1200.2 


13.785 ±0.031 






N II 


1083.9 


> 15.234 






Si II 






14.179 ±0.012 


-1.741 ±0.101 


Si III 


1206.5 


> 13.545 






S III 


1012.5 


14.562 ± 0.043 






S III 


1190.2 


13.588 ± 0.089 






Ar I 


1066.6 


< 13.256 


< 13.256 


< -1.624 


Fe II 






14.030 ± 0.032 


-1.830 ±0.105 


Fe III 


1122.5 


< 13.229 
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Fig. 10. — Velocity plot of the metal-line transitions for the damped Lyo system at z = 3 toward QXO0001. 
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2.12. Q23U+12, z = 2.538 

L98 presented a measurement of N° based on their HIRES 
observations. Our data covers the same transitions as 
well as several additional profiles (Figure 9). Unfortu- 
nately, the NI 2000a transition is crowded between two 
Lya clouds and it is difficult to accurately assess the N° 
column density. Our adopted value, 7V(N°) = io 13 - 779±0 03 cm' 
is k, 0.2 dex larger than the value reported by L98. This 
difference is substantial and is probably related to mea- 
surements of the continuum level near the NI 1200a pro- 
file. Unlike L98, our observations include the NI 1334 
triplet and we find an N° column density from these lines 
in good agreement with the NI 1200a transition. Ta- 
ble 10 lists all of the column density measurements re- 
lated to our N/a analysis. This includes a reported upper 
limit to iV(Fe ++ ) where the observed profile is most likely 
blended with several coincident absorption lines. This con- 
clusion stems from the observed differences between the 
Felll 1122 profile and the low-ion transitions, particularly 
at v w +5 km s _1 . 



Table 12 

IONIC COLUMN DENSITIES: QXO0001, z = 3 



Table 11 

IONIC COLUMN DENSITIES: Q2348-14, 



2.279 



Ion 


A 


AODM 


-^adopt 


[X/Hj 


HI 


1215.7 


20.560 ± 0.075 






N I 


1200.2 


< 13.235 


< 13.235 


< -3.295 


Al II 






12.657 ± 0.006 


-2.393 ±0.075 


Al III 


1854.7 


12.726 ± 0.012 






Al III 


1862.7 


12.613 ± 0.027 






Si II 






14.203 ±0.020 


-1.917 ±0.078 


Si III 


1206.5 


> 13.499 






S II 






13.725 ±0.119 


-2.035 ±0.141 


S III 


1190.2 


< 13.945 







-14, z = 2.279 



2.13. 

Most of the important transitions for this damped Lya 
system were presented in Prochaska & Wolfe (1999) in- 
cluding the NI 1200a transition on which our 7V(N°) limit 
is based. Table 11 lists the column densities taken from 
that paper with revised oscillator strengths where appli- 
cable. We also report an upper limit to 7V(S ++ ) based on 
the partially blended S III 1190 profile. 

2.14. QXO0001, z = 3 

Figure 10 presents metal-line profiles for the damped 
Lya system at z = 3 toward QXO0001 (logiV(HI) = 
20.7±0.05) which we will discuss further in a future paper. 
Our observations include several intermediate-ion transi- 
tions which we expect are blended with coincident Lya 
clouds. Table 12 provides the ionic column density mea- 
surements and limits for all of these transitions. Because 
the Si III 1206 transition arises in one wing of the damped 
Lya profile, we first renormalizcd the continuum surround- 
ing the Si III 1206 profile to account for this depression. 
The profile is slightly saturated and one may consider our 
7V(Si ++ ) value to be formally a lower limit. 



Ion 


A 


AODM 


N A t 

ly adopt 


IX/HI 


HI 


1215.7 


20 700 ± 050 






C III 


977.0 


> 14.171 






N I 


1134.1 


< 14.127 


1 3 31 5 -I - 040 


Q ore _l_ Q 064 


2 N I 


1200.2 


1 3 325 ± 043 






' N I 


1200.7 


1 3 272 ± 095 






N II 


1083.9 


< 13.637 






N III 


989.7 


> 14.738 






O I 


976.4 


15.766 ±0.024 


15.766 ±0.024 


-1.804 ±0.055 


O I 


1039.2 


> 15.745 






Si II 


989.8 


> 14.685 


14.452 ± 0.006 


-1.808 ±0.050 


Si II 


1193.2 


> 14.280 






Si II 


1304.3 


14.453 ± 0.006 






Si III 


1206.5 


13.47 ±0.10 


13.47 ±0.10 




S III 


1012.5 


< 13.460 






Ar I 


1048.2 


< 13.384 


< 13.384 


< -1.836 


Fe II 


1063.9 


< 15.092 


< 15.092 


< -1.108 


Fc III 


1122.5 


< 12.965 







3. DETERMINING [N/a] AND [a/H] 

The previous section presented all of the systems with 
observed NI transitions in our HIRES database. For the 
analysis which follows, we have synthesized these obser- 
vations with all other high precision nitrogen measure- 
ments for damped Lya systems at z > 1.5. This in- 
cludes the important sample introduced by L98, the sys- 
tem toward Q0201+11 by Ellison et al. (2001), and sev- 
eral new measurements (Molaro ct al. 2000; Levshakov 
et al. 2002; Dessauges-Zavadsky, Prochaska, & D'Odorico 
2002b; Lopez et al. 2002) obtained with the UVES-VLT 
spectrograph (Dekker et al. 2000). Tabic 13 summarizes 
all of these observations. Column 6 of the table indicates 
the expected importance of ionization corrections as de- 
scribed below and column 7 indicates whether the system 
is included in the analysis of § 4. In general, those sys- 
tems where we expect ionization corrections will be sig- 
nificant but could not reasonably assess their value were 
eliminated. 

Given the appropriate wavelength coverage, high reso- 
lution observations of the damped systems can generally 
provide precise and accurate ionic column densities for 
both nitrogen and an a-element. Before converting these 
gas-phase measurements into elemental abundances and 
N/a ratios, however, one must consider several factors: 
(1) identifying the 'best' a-element; (2) dust depletion; 
and (3) photoionization. At present, we believe the first 
two issues have relatively minor impact on derivations of 
N/a and a/H in the damped systems. In previous analy- 
ses and this work too, authors have relied on Si to serve as 
the a-element because it is the most readily observed. This 
raises some concern because Si is expected to be synthe- 
sized in less massive progenitors of Type II SN than O or 
Mg (Woosley & Weaver 1995) and may even be substan- 
tially produced in Type la SN (Matteucci & Recchi 2001). 
Therefore, reporting [N/a] as [N/Si] may complicate com- 
parisons against theoretical models which generally reflect 
nucleosynthesis in the most massive stars. Empirically, 
however, we and others have found that [O/Si] ~ and 
[S/Si] « in the few damped systems where Si and O or 
S have been determined. Therefore, we believe the differ- 
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Table 13 
SUMMARY 



QSO 


z abs 


AT(HI) 


[a/Hj* 


[N/aj 


IC C 


Included 


Rcf 


Q0000-2619 


3.390 


21.41 


-1.91 ±0.08 


-0.74 ± 0.04 


Y 


Y 


1,2,3 


PH957 


2.309 


21.37 


-1.46 6 ±0.08 


-0.85 ± 0.02 


Y 


Y 


4,5 


Q0201+11 


3.387 


21.26 


-1.25 6 ±0.15 


-0.65 ± 0.16 


? 


Y 


6 


Q0201+36 


2.463 


20.38 


-0.41 ±0.05 


> -0.94 


Y 


Y 


7,8 


J0307-4945 


4.468 


20.67 


-1.55 ±0.12 




N 


N 


9 


Q0336-01 


3.062 


21.20 


-1.41 6 ±0.10 


> -0.73 


Y 


Y 


8 


Q0347-38 


3.025 


20.63 


-1.17 ±0.03 


-0.54 ± 0.04 


Y 


Y 


8, 10 


HS0741+4741 


3.017 


20.48 


-1.69 ±0.10 


-0.79 ±0.01 


Y 


Y 


8 


Q0930+28 


3.235 


20.30 


-1.97 ±0.10 


-0.56 ±0.02 


N 


Y 


8, 11 


Q1055+46 


3.317 


20.34 


-1.65 ±0.15 


< -0.57 


? 


Y 


11 


BR1202-07 


4.383 


20.60 


-1.81 ±0.14 


< -0.46 


? 


Y 


11 


Q1223+17 


2.466 


21.50 


-1.59 ±0.10 


-1.05 ± 0.18 


N 


N 


8 


Q1331+17 


1.776 


21.18 


-1.45 ±0.04 


-0.49 ±0.11 


? 


Y 


5 


Q1425+6039 


2.827 


20.30 


> -1.03 


< -0.53 


N 


N 


8, 11 


Q1759+75 


2.625 


20.76 


-0.89 ±0.01 


-0.68 ± 0.02 


N 


Y 


12 


Q1946+7658 


2.844 


20.27 


-2.23 ±0.06 


-1.43 ± 0.04 


Y 


Y 


8, 11 


Q2212-1626 


3.662 


20.20 


-1.90 ±0.08 


< -0.69 


? 


Y 


11 


BR2237-0607 


4.080 


20.52 


-1.87 ±0.11 


< -0.34 


? 


Y 


11 


HE2243-6031 


2.330 


20.67 


-0.87 ±0.03 


-0.89 ± 0.02 


? 


Y 


13 


Q2343+12 


2.431 


20.34 


-0.54 ±0.10 


-1.10 ± 0.09 


N 


N 


5 


Q2344+12 


2.538 


20.36 


-1.74 ±0.10 


-0.81 ±0.03 


Y 


Y 


8 


Q2348-14 


2.279 


20.56 


-1.92 ±0.08 


< -1.38 


N 


Y 


2,8,14 


QXO0001 


3 


20.70 


-1.81 ±0.05 


-1.45 ± 0.04 


Y 


Y 


8 



"Assumes Si unless noted 
b Sulfur 

c See below for a detailed discussion of ionization corrections. 
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ence between O and Si leads to an uncertainty in [N/a] of 
less than 0.1 dex. 

We expect that dust depletion has an even smaller effect 
on the [N/a] and [a/H] values than uncertainties between 
the various a-elements. Although Si is significantly de- 
pleted along many ISM sightlines (e.g. Savage & Sembach 
1996), the depletion levels in the damped Lya systems (as 
inferred from Zn/Fe or S/Fe) are too small for the dif- 
ferential depletion of Si/N or Si/H to be important (e.g. 
Vladilo 1998). This is particularly true for the systems 
presented here where the Si/Zn and Si/S ratios are very 
nearly solar (Prochaska & Wolfe 2002; but see Ledoux et 
al. 2002 for a counter-example). The largest corrections 
one could derive for Si relative to N and H is +0.2 dex 
which would imply slightly lower N/a ratios and a higher 
a/H value. 

The major sources of systematic error in converting N° , 
Si + , H , and O column densities into elemental abun- 
dances are ionization corrections. Consider first the ion- 
ization balance of O , 0+, N°, N+, H°, and H+. Charge- 
exchange reactions between (O , H°) and (N°, H°) help 
to regulate the ionization fractions of these ions such that 
H°/H « N°/N « O /O when the volume density is large 
and the ionizing flux is modest. If the ratio of ionizing 
flux to volume density - the ionization parameter U - 
is large enough, then this set of equalities breaks down. 
Both O and N° have larger cross-sections to photons with 
hv > 2 Ryd than H° (Sofia & Jenkins 1998) and these ions 
will be over-ionized: N°/N < H°/H and O /O < H°/H. 
The likelihood of over-ionization for Ar is even greater as 
found in the Galactic ISM (Sofia & Jenkins 1998) and in 
damped systems (P02). In contrast, Si + can be present in 
highly ionized regions such that Si + /Si > H°/H for nearly 
all physical conditions. The net result is that Si + /H° will 
overestimate Si/H and N°/Si + will underestimate N/Si if 
ionization is significant in the damped system, i.e., when 
\ogU ;» — 3 (see below). The magnitude of these correc- 
tions ranges from a few hundredths dex in a very neutral 
gas to several tenths dex in an ionized gas. To accurately 
assess these ionization corrections, one must examine pho- 
toionization diagnostics within each damped system. 

Historically, researchers have assumed that ionization 
corrections are negligible for the damped systems. This 
presumption is physically motivated by the very large op- 
tical depth at energies hv > 1 Ryd implied by a sight- 
line with N(Rl) > 2 x 10 20 cm-2 (Viegas 1994). De- 
tailed theoretical models have drawn similar conclusions 
for the majority of elements observed in the damped sys- 
tems (Howk & Sembach 1999; Vladilo et al. 2001). More 
recently, however, Prochaska et al. (2002) presented evi- 
dence for significant photoionization in w 50% of the gas 
associated with the damped Lya system at z = 2.62 to- 
ward GB1759+75. For this gas, the ionization corrections 
are important for ions like N° and Si + relative to H°. The 
correction for N°/Si + in this partially ionized gas is large: 
[N/a] = [N°/Si+] + 0.3 dex. Although the authors stress 
that this particular damped system has several charac- 
teristics which separate it from the majority of damped 
systems, its properties highlight the importance of assess- 
ing the ionization state of each damped system in order to 
convert the observed N°/Si+, N°/S+, 0°/H°, S+/H , and 
Si+/H° ratios to N/a and a/H. 



In Figure 11, we present a series of ionic ratios which 
qualitatively assess the level of ionization corrections for 
gas with a range of ionization states. These values were 
calculated with the CLOUDY software package (v95b; Fer- 
land 2001) assuming a Haardt-Madau extragalactic back- 
ground radiation field at z = 2.5 (Haardt & Madau 1996). 
In panel (a) we plot N°/Si + and Si + /H° versus the ioniza- 
tion parameter, U = </>/ (cuh) where <p is the surface flux of 
ionizing photons with hv > 1 Ryd. The vertical dotted line 
at log U — —3.5 demarcates the transition from corrections 
of less than 0.2 dex for [N°/Si + ] to greater than 0.2 dex. 
At this U value, panel (b) shows: [Fe++/Fe+] w -0.8, 
[Si++/Si+] » -1.1, [S++/S+] w -0.35, [Ar°/Si+] ss -0.7, 
[A1++/A1+] w -0.6, and [N+/N ] w -0.5. Systems which 
exhibit significantly lower values for these ratios (higher 
values for Ar°/Si + ) should be predominantly neutral and 
require small corrections for [N/a] and [a/H]. Similarly, 
systems with larger values of these ratios are presumably 
partially ionized and may require significant ionization cor- 
rections. In the following sub-sections, we will describe the 
ionization diagnostics for each of the damped systems in 
the full nitrogen sample and evaluate the magnitude of 
ionization corrections. 

3.1. Q0000-26, z = 3.390 

Molaro et al. (2001) presented Ar° measurements for 
this damped Lya system which indicate a nearly solar 
Ar/Si ratio. We agree with their argument that the gas in 
this damped Lya system is therefore predominantly neu- 
tral. 

3.2. PH957, z = 2.309 

Dessauges-Zavadsky, Prochaska, & D'Odorico (2002b) 
will present a detailed analysis of the photoionization state 
of this damped system. Their assessment, based on [Ar°/S + ] ~ 
—0.2 and observations of the Felll 1122 transition, is that 
ionization corrections for N° and S + are small. 

3.3. Q0201+11, z = 3.387 

There is no photoionization diagnostic for this damped 
system, but we expect that its HI column density is large 
enough that photoionization corrections should be small. 

3.4. Q0201+36, z = 2463 

Our new data on this system covers the Felll 1122 tran- 
sition and provides a measurement of Ar°. The [Ar°/Si + ] 
value of —0.42 ± 0.1 dex suggests the possibility of mild 
photoionization while the observed Fe ++ /Fe + ratio of — 1 dex 
indicates such effects are likely to be small. Altogether we 
expect ionization corrections of 0.1 dex or less and assume 
no correction. For this damped system, we measure ions 
for two a-elements, Si + and S + whose abundances rela- 
tive to solar are in reasonably good agreement. For [a/H], 
we adopt the [Si/H] value because the Sill 1808 profile 
is outside the Lya forest and therefore less likely to be 
contaminated by line-blending. 

3.5. J0307-4945, z = 4466 

Dessauges-Zavadsky et al. (2001) presented a detailed 
analysis of this damped system from their UVES obser- 
vations which included measurements of N, Si, O and Fe 
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Fig. 11. — Ionic ratios predicted for a plane-parallel slab of gas with JV(HI) = 10 20 cm — 2 for a range of ionization parameter U. These 
calculations were carried out with the CLOUDY software package. The dotted vertical line indicates the transition from a predominantly 
neutral gas where ionization corrections to [N°/Si+] are small to a partially ionized gas where ionization corrections are significant. Provided 
the appropriate observations, the calculations presented in this figure can be used to roughly assess the ionization corrections for a given 
damped Ly« system. 



among other elements. The authors presented N° column 
density measurements for the strongest velocity compo- 
nents of this kinematically complex system and discussed 
the relatively low N/a value. We are concerned, however, 
that this gas may require significant ionization corrections. 
We note that for those components where the O I profiles 
were unsaturated, the measured O /Si + ratios are signif- 
icantly sub-solar, in some cases by over 1 dex. This is 
a clear signature of photoionization in these weaker com- 
ponents where NI is undetected. Unfortunately, the OI 
profiles are heavily saturated in those components where 
N I was measured and it is impossible to accurately assess 
the photoionization state of the gas from the data pre- 
sented in Dessauges-Zavadsky et al. (2001). Although we 
consider it a reasonable possibility that this system does 
exhibit a low N/a value (consistent with the low nitro- 
gen sub-sample discussed in § 4) , we choose to eliminate it 
from our sample on the grounds that photoionization may 
be important. 



3.6. Q0336-01, z =3.062 

Our observations of the damped Lya system toward 
Q0336-01 provide a lower limit on the Ar°/S + ratio and 
place upper limits on the N + /N° and Fe ++ /Fe + . The 
nearly solar [Ar°/S + ] ratio and the relatively small N + /N° 
and Fe ++ /Fe + limits all indicate this system is predomi- 
nantly neutral. 

3.7. Q0347-38, z = 3.025 

To assess the ionization state of this damped system we 
focus on the Ar° measurements made by Levshakov et al. 
(2002). Their analysis indicates [Ar°/Si+] > -0.2 imply- 
ing this system is neutral. This conclusion is supported by 
our measured ratio of Fc ++ /Fc + = —1.4 dex which may 
be considered an upper limit (§ 2.4). We also report a 
S ++ /S + value of < —0.8 dex which is consistent with this 
conclusion. 

It is worth noting, that the C III 977, Felll 1122, S III 1012, 
and Si III 1206 profiles all exhibit significant absorption at 
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v > 20 km s _1 . Although this absorption may result from 
coincident Lya clouds, we contend the absorption is more 
likely due to a significantly ionized gas associated with the 
damped system. This component does not give rise to sig- 
nificant low-ion absorption, however, and therefore has no 
effect on our N/a analysis. 

3.8. HS0741+47, z = 3.017 

Our spectrum of HS0741+47 includes coverage of several 
ionization diagnostics for this damped Lya system includ- 
ing Arl, Felll, and Nil transitions. Unfortunately, the 
Felll 1122 and Nil 1083 transitions appear to be blended 
with coincident Lya forest clouds (Figure 4) and the re- 
sulting upper limits are not too meaningful. On the other 
hand, [Ar°/Si + ] = —0.18 and we contend the system is 
primarily neutral with ionization corrections for [N/a] and 
[a/H] less than 0.1 dex. 

3.9. Q0930+28, z = 3.235 

The strongest constraints on the ionization state of this 
system are the N+/N and Si++/Si+ ratios. The Nil 1083 
profile may be mildly blended with a coincident Lya for- 
est cloud, but we expect N + /N° w —0.2 dex which im- 
plies a modest ionization correction for nitrogen. Simi- 
larly, the saturated Si III 1206 profile yields a lower limit 
to Si ++ /Si + of —0.25 dex consistent with this system be- 
ing partially ionized. Altogether, we contend the N/a and 
a/H values would require significant ionization corrections 
and we do not include this system in our analysis. In 
passing, however, we note that the uncorrected [N/a] and 
[a/H] values are typical for the majority of our sample. 

3.10. Q1055+46, z = 3.317 

We have no ionization diagnostics for this system. 

3.11. BR1202-07, z = 4.383 

We have no ionization diagnostics for this system. 

3.12. Q1223+17, z = 2.466 

As described in § 2.7, we have been unable to measure an 
accurate N° column density from our dataset and have not 
included this system. In passing, we note the Felll and 
N II profiles indicate the gas is predominantly neutral at 
v k, km s _1 but suggest the gas at v w —70 km s _1 could 
be partially ionized. The latter point hinges on whether 
the feature observed at v « —70 km s _1 in the Felll 1122 
profile is the result of a blend. 

3.13. Q1331+17, z = 1.776 

A complete analysis of this system's ionization state 
will be presented in Dessauges-Zavadsky, Prochaska, & 
D'Odorico (2002b). Their UVES observations of Fe III 1122 
indicate the gas giving rise to the N I profiles is predomi- 
nantly neutral. 

3.14. Q 1425+ 60, z = 2.827 

This system exhibits several diagnostics which indicate 
photoionization is likely to be important. In particular, we 
find [Ar°/Si+] < -0.5 dex and N+/N > dex. Although 



the Nil 1083 profile is mildly saturated and therefore re- 
sembles a typical Lya cloud, we are confident it is not a 
coincident absorption line because the profile also exhibits 
a velocity component at v w —90 km s _1 which coincides 
with a similar feature in the low-ion profiles. Analysing 
this single velocity component, we find N + /N° w dex 
consistent with the total N(N+) and AT(N°) values. Alto- 
gether we expect the measured N° column density signif- 
icantly underestimates the nitrogen elemental abundance. 
We choose to remove it from the analysis noting, however, 
that it does not appear to exhibit a low N/a ratio. In fact, 
if ionization corrections to N° and Si + are significant, this 
system may exhibit a nearly solar N/a ratio. 

3.15. GB1759+ 75, z = 2.625 

A detailed photoionization treatment of this damped 
system has been presented by P02. We concluded that 
roughly half of the gas is partially ionized with N°/Si + re- 
quiring significant ionization corrections. The remaining 
gas, meanwhile, is predominantly neutral. We are able to 
estimate the ionization corrections to [N/a] and [a/H] for 
the entire system and have therefore included it within our 
analysis. The corrections we adopt are +0.15 dex to [N/a] 
and —0.1 dex to [a/H] as measured from [Si/H]. 

3.16. Q1946+ 76, z = 2.844 

This well studied damped Lya system has a relatively 
low HI column density and one of the lowest N°/Si + ratios 
observed. These factors raise concerns that its gas may be 
primarily ionized and that the low inferred N/a value has 
resulted from photoionization. This concern is supported 
by the slightly sub-solar O /Si + ratio, but we place a strict 
upper limit on Si ++ /Si + < —0.75 dex which suggests ion- 
ization corrections are not extreme. Unfortunately, the 
Nil 1083 profile is blended and does not provide a useful 
limit on N + /N°. Nevertheless, the system may require an 
ionization correction for [N/a] (see Figure 11) and we have 
remaining concerns on its ionization state. In the follow- 
ing, however, we proceed without adopting any correction. 

3.17. Q2212-16, z = 3.662 

There is no ionization diagnostic for this damped sys- 
tem. 

3.18. BR2237-06, z = 4.O8O 

There is no ionization diagnostic for this damped sys- 
tem. 

3.19. HE2243- 60, z = 3.149 

This new damped Lya system (Lopez et al. 2002) ex- 
hibits competing photoionization diagnostics. Although 
the system shows a significantly sub-solar Ar°/Si + ratio, 
[Ar°/Si+] = -0.5, Lopez et al. (2002) examined the A1+ 
and Fe + column densities and argued that ionization cor- 
rections are likely to be small. We proceed under the as- 
sumption that ionization corrections to N/a are small. 

3.20. Q2343+12, z = 2.431 

Dessauges-Zavadsky, Prochaska, & D'Odorico (2002b) 
will present a full analysis on the photoionization of this 
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system. We note here that they find a very low Ar°/Si + 
ratio and strong N + absorption indicating photoionization 
corrections will be important. As such, we have not in- 
cluded this system in our analysis. 

3.21. Q23U+12, z = 2.538 

Although our observations cover several ionization diag- 
nostics, the only valuable ratio is Fc ++ /Fc + ; Nil 1083 is 
severely blended and the Ar° transitions provide only an 
upper limit to 7V(Ar°). The upper limit on Fe ++ /Fe + < 
—0.8 suggests the system is probably neutral but the N° 
and Si + column densities could be subject to small (< 
0.1 dcx) ionization corrections. In the following, we adopt 
no corrections. 

3.22. Q2348-U, z = 2.279 

The very low upper limit on N°/Si + for this damped 
system raises concerns that photoionization may be im- 
portant. Unfortunately, the existing spectra do not cover 
the FcIII or Nil profiles. We do measure Al++/Al+ = 
—0.55 dex but this diagnostic is relatively insensitive to the 
photoionization state (Figure 11) and does not provide a 
particularly meaningful constraint. We also place a lower 
limit on Si ++ /Si + of —0.7 dex which may be influenced 
by a line-blend in the Lya forest. Unfortunately, these 
diagnostics are inconclusive. Most worrisome is the fact 
that the Si IV and C IV profiles exhibit significant features 
at essentially identical velocity as the low-ion absorption 
indicating ionized gas is at this velocity. As discussed in 
P02 and Wolfe & Prochaska (2000), it is unlikely that the 
gas responsible for the high-ion profiles is associated with 
the low-ion profiles in damped Lya systems. Nevertheless, 
we have remaining concerns regarding the ionization state 
of this system. On the other hand, the N° column density 
is an upper limit and we are reasonably confident that the 
N/a is very low in this damped system. 

3.23. QXO0001, z = 3 

Although our observations provide several ionization di- 
agnostics for this system, it is difficult to draw a definite 
conclusion on its ionization state. The N II 1083 profile is 
blended and places a large upper limit on N + /N°. The 
Ar I transition is also blended and provides an upper limit 
[Ar°/Si + ] < dex which is consistent with both a neu- 
tral and ionized gas. Meanwhile, the partially saturated 
Si III 1206 transition implies Si ++ /Si + > — 1 and we place 
an upper limit on Fc ++ /Fc + of —1 dex. These latter 
two diagnostics suggest the system is primarily neutral, 
but a mild ionization correction to N/a is allowed (up to 
0.2 dex). We proceed assuming a +0.1 dex correction to 
[N/a] and no correction to [a/H]. 

4. ANALYSIS AND DISCUSSION 

4.1. What The Data Reveal 

Abundance results from Table 13 are presented graph- 
ically in Figure 12. This figure shows the abundance ra- 
tios [N/a] versus [a/H] for a large sample of H II regions 
and blue compact galaxies compiled by Henry & Worthey 
(see their Figure 6 along with the caption, which identifies 
the numerous data sources). To their sample we have now 



added our DLA results, which are shown in the figure with 
open circles, triangles (limits), and error bars. 

The quantity [a/H] in the figure is a gauge of metallic- 
ity and refers to those heavy elements which are produced 
through subsequent additions of alpha nuclei. In the case 
of the Henry & Worthey sample, a corresponds to oxygen, 
while for the objects in our sample, it refers to silicon in 
most cases, and sulfur in a few cases (Table 13). Metallic- 
ity is in turn a measure of the extent to which the observed 
gas, which presumably began as pristine and metal-poor 
in these objects, has been processed through stars, since 
metal abundances rise as nuclear products from evolved 
stars are added to the gas. Metallicity is directly related 
to the total amount of star formation occurring since the 
system formed, and therefore values increase with time 
with a slope determined by the star formation rate and 
are tempered by the accretion of primordial gas from the 
surrounding IGM. On the other hand, [N/a] measures the 
differential change in nitrogen relative to a-elements as 
processing goes forward. It is sensitive to the stellar yields 
for each element as well as the form of the initial mass 
function, but relatively insensitive to the detailed history 
of the star formation rate. 

As we study the global pattern created by the emission 
line objects in Figure 12, we identify a lower-bound en- 
velope behind which points are concentrated in number 
densities that fall off with increasing distance from the en- 
velope. Referring to the discussion in § 1, the plateau 
in the envelope at low metallicity that is displayed by 
the metal-poor galaxies (Izotov & Thuan 1999; Kobul- 
nicky & Skillman 1996) is explained by the dominance 
of primary (metal-insensitive) nitrogen production at low 
metallicity. Likewise, the upturn at higher metallicities 
can be attributed to increasing contributions from sec- 
ondary (metal-sensitive) production. Currently, the causes 
of the scatter behind the envelope are not well-understood, 
although several theories have been proposed (Garnett 
1990; Pilyugin 1993). 

Most of the data points in Figure 12 collectively repre- 
sent emission line regions in many different host galaxies 
of several different galaxy morphologies. Therefore, the 
fact that we see a pattern formed by this diverse group 
of objects suggests that, to first order, universal synthe- 
sis patterns of stars, as communicated to the environment 
through their IMF- weighted yields, are more important 
than local factors such as the star formation rate, infall, 
outflow, etc. for determining how nitrogen and a-element 
abundances evolve with respect to each other. 

Now consider the damped Lya systems plotted in Fig- 
ure 12. The majority follow a lower metallicity extension 
of the N/a plateau expressed by the emission line regions, 
and unlike previous studies (e.g. L98), we observe minimal 
scatter in this group. The DLA 'plateau' values exhibit an 
rms scatter < 0.2 dcx which is remarkably small given the 
systematic uncertainties described in § 3. Furthermore, 
all of the limits in this sub-sample are consistent with an 
[N/a] value of —0.75 dex. In contrast with the emission 
line regions, we identify three low nitrogen damped sys- 
tems (LN-DLA) which represent a region of the plot that 
is isolated from and discontinuous with the main pattern. 
Similar to the major DLA population, these LN-DLA also 
exhibit minimal dispersion in their N/a values. Therefore, 
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Fig. 12. — [N/a] versus [a/H] for a large sample of H II regions and stars compiled by Henry, Edmunds, & Koppcn (2000) and Henry 
& Worthey (1999) and the damped Lyman alpha systems studied in this paper. Symbols for all but the DLAs are defined in Henry et al. 
(2000). The DLAs are shown with open circles and error bars. Up triangles represent lower limits, down triangles upper limits. 



the DLA systems form a bimodal distribution where the 
majority lie on the metal-poor plateau and a sub-sample 
exhibit significantly lower N/a values. 

This bimodality of the DLA N/a measurements was not 
apparent in previous studies. In part, the difference is 
our reliance on higher quality observations, but of simi- 
lar importance is our greater appreciation of the effects 
of photoionization. Our current sample is small enough, 
however, that we do not have complete confidence in a 
bimodal distribution. It is worth noting, however, that 
the new measurements presented in Pettini et al. (2002) 
lend further support to this description. In the follow- 
ing, we will investigate the implications of the LN-DLA 
sub-sample and the impact of a bimodal distribution on 
various models of N production and star formation. 

4.2. Model Interpretations 

Several studies cited in §1 indicate that IMS produce 
the bulk of the nitrogen in the universe. Therefore, we 
expect that the LN-DLA arise as the result of a reduced 
IMS contribution to nitrogen buildup. Here, we list four 
ways that such a reduction could occur. 

1. The LN-DLAs are observed during a transient 
period occurring within 250 Myr of a star burst 
and before the slower evolving IMS have released 
their nitrogen into the interstellar medium of these 
objects. 



2. Nitrogen yields of IMS significantly decline with a 
decrease in metallicity. 

3. The observed [N/a] reflects the nucleosynthesis 
characteristics of Population III stars. 

4. The IMF in these objects is either truncated at the 
low mass end or the slope is significantly flatter 
than the commonly observed form, resulting in 

a much larger ratio of massive stars to IMS and 
reducing the contribution of IMS to nitrogen levels. 

Because the first three of these have been put forth by 
other authors, we only describe them briefly and then turn 
our attention to our fourth point. For each option, we 
examine the scenario under the constraint that it naturally 
reproduces a bimodal distribution of N/a values. We will 
find that this constraint severely challenges the first three 
models and, therefore, our results favor a top heavy IMF 
mode of N production in the LN-DLA. 

The first hypothesis has been considered previously (L98; 
HEK00; Pettini et al. 2002) and is a reasonable possibility 
when one considers that at z ~ 3, the age of the interven- 
ing galaxies must be less than «2 Gyr, i.e., the age of the 
universe. It is possible that we are observing the damped 
Lya systems during very early formation times when their 
IMS would not have fully expelled their N products. This 
was shown to be a viable explanation by HEK00, who esti- 
mated the time delay in IMS nitrogen release to be about 
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Fig. 13. — [N/q] (solid curve) and [a/H] (dotted curve) versus 
time in Myr for model B of HEK00. The horizontal dashed lines 
indicate ccntroid values for observed [N/a] and [a/H] associated 
with the LN-DLA objects displayed in Figure 12. 



250 Myr 5 . The challenge for this scenario is revealed by 
Figure 13 which plots [N/a] and [a/H] versus time in Myr 
for model B of HEK00, the model which passes nearest the 
location of the LN-DLAs. Ovcrplotted on the theoretical 
curves are the N/a and a/H values (horizontal dashed 
lines) of the LN-DLA sample. The plots show that the 
required levels are not reached at the same time: [N/a] 
reaches a level of —1.4 at the same time that [a/H] is 
about 0.5 below —2.0 dex. This inconsistency could be 
overcome, however, by fine-tuning the star formation effi- 
ciency More important is that model B and its variants 
predict the gas will spend as much time at [N/a] ~ —1.4 as 
it will at —1.2, —1.6 dex, and other low values. Therefore, 
one would expect the LN-DLA to uniformly fill the region 
below the higher N/a plateau expressed by the majority 
of DLA. The identification of the three LN-DLA at nearly 
the same N/a value must be explained as a coincidence, 
i.e., this scenario cannot naturally reproduce a bimodal 
distribution. 

In considering option 2 above, we employ a simple closed- 
box model of chemical evolution in order to estimate the 
ratio of nitrogen to oxygen yields needed to explain the 
observed [N/a] of -1.4 of the LN-DLA. Recall that the 
metallicity in a simple model can be expressed as Z x = 
—yxlnfi, where Z x is the mass fraction of element x, y x 
is the total yield of a population of stars characterized by 
a Salpeter IMF per unit of stellar material forever locked 
up in remnants, and \i is the gas fraction (Tinsley 1980). 
From the above observed value of [N/a] and the solar 
value for N/O of -0.9 (Grevesse & Noels 1996) we infer 
that Z^jZo = Vn/vo = 0.005. HEK00 integrated several 
sets of published yields over a Salpeter IMF, and we can 
now compare their results for low metallicity directly with 

5 Siess, Livio, & Lattanzio (2002) have recently calculated IMS evo- 
lution models in which they find that the stellar mass threshold for N 
production drops at lower metallicity, thus increasing the timescale 
for its production by IMS even further. 



this observed value. Such a comparison shows that current 
published yield predictions of N/O are about two orders of 
magnitude above the level implied by the LN-DLA obser- 
vations. Their results appear to rule out option 2 above, 
although published IMS yield calculations so far do not 
extend down to metallicities below 0.001. Furthermore, 
there is no indication from published IMS nitrogen yields 
that they are highly sensitive to metallicity (van den Hoek 
& Groenewegen 1997; Marigo 2001). Finally, it would be 
very difficult to explain the bimodality of [N/a] values at 
the same [a/H] metallicity in this scenario. 

For option 3 we compare in a similar way the same 
simple model results to predicted yield ratios. Heger & 
Woosley (2002) have calculated models of massive Pop III 
stars ranging in mass from 140 to 260 M Q and predicted 
yields for 14 N and 16 O. Simply by inspecting their results 
we can see that values of N/O resulting from their mod- 
els range from about 10 -6 in the least massive stars to 
10 -7 in the most massive ones. In addition, Woosley & 
Weaver (1995) predict N and O production by zero metal- 
licity stars between 12-40 M and find values of N/O of 
a few times 10~ 4 , where again we have used the HEK00 
results for yields integrated over a Salpeter IMF for com- 
parison. These values are lower than necessary to ex- 
plain LN-DLAs. However, in their Population III stel- 
lar models Umeda, Nomoto, & Nakamura (2000) predict 
14 N/ 24 Mg values of roughly 1/4 solar, which they say im- 
plies a large contribution of Population III stars to the 
early nitrogen buildup in their yield ratios which are con- 
sistent with these. In summary, while the relevance of 
Population III stars is presently unclear, these stars could 
ultimately prove important for the issue of early nitro- 
gen production. Furthermore, if one could introduce a 
Pop III model which reproduced the observed LN-DLA 
population, it could reasonably account for the bimodal 
distribution where the higher N/a DLA represent a sec- 
ond generation of star formation. 

For the remainder of the discussion, we explore the fourth 
hypothesis: star formation is either characterized by (1) a 
top-heavy initial mass function at early times and/or low 
metallicities, such that there is a relative deficiency of IMS; 
or (2) a stellar population that is truncated below a certain 
mass threshold. In either case, the observed [N/a] in the 
LN-DLAs is dominated by the yield patterns of massive 
stars. 

To test this hypothesis, we employ the numerical code 
of HEK00 and calculate two types of chemical evolution 
models, both of which assume a closed system, i.e. no 
matter exchange with the environment. In the first set of 
models we ignore evolutionary time lags among stars of dif- 
ferent progenitor masses (instantaneous recycling approxi- 
mation) , while in the second set we relax this assumption. 
Otherwise the model types are identical. The numerical 
code we employed is described in detail in HEK00. The 
only alterations we made relate to the star formation rate 
(eq. 18 in HEK00), where we used 1.5 rather than 2.0 for 
the exponent and set the efficiency factor to a constant 
value of 0.003 throughout. 

We tested two sets of yields comprising the following 
combinations: (1) Woosley & Weaver (1995; massive stars) 
and van den Hoek & Groenewegen (1997; IMS); and (2) Porti- 
nari et al. (1998; massive stars) and Marigo (2001; IMS). 
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Fig. 14. — (a) Model results for set 1 stellar yields displayed with data points of Figure 12. Light solid curves represent results for the two 
different IMF slopes indicated and when instantaneous recycling is assumed (IRA+). Heavy solid lines show tracks for same models except 
that the instantaneous recycling assumption is relaxed (IRA-). The dashed line shows the model result when the IMF is truncated below 
7.50 Mq and instantaneous recycling is assumed, (b) Same as (a) but for set 2 stellar yields. 



Our calculations, therefore, accounted for the effects of 
progenitor mass and metallicity on the yields. The IMF 
was taken to be a power law form, the index of which we 
varied so as to alter the relative mix of massive stars and 
IMS. For the models which assumed instantaneous recy- 
cling, we calculated an additional model subset for each 
yield set using a Salpeter IMF (index=1.35) and exper- 
imented with different lower mass threshold values until 
finding one which satisfactorily reproduced the observed 
[N/a] value. Besides the choice of yield set, the only pa- 
rameters allowed to change in our models were the slope 
and lower mass threshold of the IMF. 

The model results assuming instantaneous recycling are 
shown with solid lines in Figures 14a, b, where each graph 
represents models associated with one yield set (see figure 
legends for line- type definitions) . Lines are labeled accord- 
ing to IMF slope or cutoff mass 6 . Because of the instan- 
taneous recycling assumption, no time lag exists between 
IMS and massive star release of N into the interstellar 
medium, and thus the tracks for these models are horizon- 
tal lines. In Figure 14a, where we employed the set (1) 
yields, the positions of the three LN-DLAs are consistent 
with a system whose stellar population possessed a top- 
heavy IMF with a slope of roughly 0.10 or with a Salpeter 
IMF but truncated below 7.5 M©. In Figure 14b, using 
set (2) yields, the LN-DLAs are again consistent with an 
IMF slope of 0.10 or a mass cutoff of 4.8 M Q . Note that 
for each yield set a slope of 1.10 produces reasonable fits to 
the plateau region for both yield sets 7 . Finally, we point 

6 Note that because we are only interested in the results at low 
metallicity, all models were stopped when [a/H]= — 1. 

7 The difference here between 1.10 and 1.35 (Salpeter slope) is not of 
immediate concern, since uncertainties in the yields themselves are 
roughly a factor of two. 



out that [N/a] is quite sensitive to the value of the mass 
cutoff in both cases, since it falls within the mass range 
where most IMS N is produced in both cases. 

Results for models without the instantaneous recycling 
approximation are shown with bold lines. The IMF slope 
in each case is identical to the model track with which the 
bold line merges close to [a/H]=-l. The gradual rise in 
[N/a] with metallicity for these models is simply due to 
the delayed contribution of IMS to the N buildup. Clearly, 
these models are still capable of reproducing DLA obser- 
vations for both the plateau and LN-DLA objects. 

Thus, our closed box models, calculated using two dif- 
ferent sets of theoretical stellar yields, appear to be consis- 
tent with the idea that LN-DLAs are systems which have 
evolved with either top-heavy or truncated IMFs. We find 
a representative IMF slope of 0.10 or a mass cutoff of 4.8 
M Q with yield set (1) and 7.5 M Q with yield set (2). This 
conclusion holds regardless of whether instantaneous recy- 
cling is assumed. The reason that the altered IMF slope or 
the truncation can explain the LN-DLAs is simply that the 
relative contribution of IMS to nitrogen evolution is sig- 
nificantly reduced, leaving the massive stars to dominate 
the production of both nitrogen and oxygen. 

The top-heavy IMF hypothesis, we believe, offers a vi- 
able alternative to ones already suggested by others and 
briefly discussed above. Furthermore, this scenario may 
also explain the bimodal distribution observed for the DLA 
measurements (see below). Indeed, the LN-DLAs may of- 
fer the first evidence that the mass spectra of some early 
stellar populations were much different than we observe 
today. 

There is currently no convincing evidence for the exis- 
tence, past or present, of an IMF which differs markedly 
from those derived from observations in local systems such 
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as the form presented by Scalo (1986). Instead, most evi- 
dence points to the IMF being invariant (e.g. Kroupa 2001) 
in space and time. However, simple arguments based upon 
a Jeans mass analysis imply that at low metallicities, star 
formation should favor higher mass objects. This factor 
alone alerts us to a potential top-heavy IMF (THIMF) at 
early times. 

Silk (1998) discusses several observations which may 
suggest a THIMF at very early epochs, although he is 
very careful to note that they are all circumstantial and 
in no way clinches the case for a THIMF. Theoretical sup- 
port for a THIMF comes from the models of Nakamura 
& Umemura (2001; 2002), who study star formation in 
extremely metal-deficient protogalactic environments and 
predict a resulting bimodal IMF with peaks around 1 and 
100 Mq. 

At present, we favor the THIMF interpretation because 
it accurately reproduces the N/a values of the LN-DLA 
and because it may be consistent with a bimodal distribu- 
tion. The latter point would hold if a primordial THIMF 
abruptly changed over to the form we observe today af- 
ter a certain threshold was passed. Bromm et al. (2001) 
have experimented with numerical simulations of the col- 
lapse of pre-enriched primordial objects and set a value of 
Z C rit = 5x 10~ 4 Z Q as the metallicity threshold when this 
changeover occurs (see also Abel, Bryan, & Norman 2000). 
In the end it seems likely that measurements of integrated 
light from early systems will prove useful in deciding be- 
tween the THIMF option and the other possibilities which 
are all consistent with a normal stellar mass spectrum. 

5. SUMMARY AND CONCLUDING REMARKS 

In this paper we presented N/a vs. a/H measurements 
of a moderate sample of damped Lya systems and identi- 
fied a possible bimodality in the N/a values at the lowest 
metallicity. While the majority of low metallicity DLA 
have N/a values consistent with low metallicity emission- 
line regions in the local universe, a sub-sample exhibits 
significantly lower N/a values (the LN-DLA). Unlike pre- 
vious works, we found minimal scatter among the mem- 
bers of each of the sub-samples and therefore argue that 
the DLA measurements form a bimodal distribution. We 
considered several scenarios which could account for these 
observations: (1) these systems are too young (< 250 Myr) 
for their intermediate-mass stars to have produced N; (2) 
a reduced production of N in intermediate-mass stars at 
low metallicity; (3) Population III nucleosynthesis; and 
(4) chemical evolution with a top-heavy IMF. Currently, 
we do not favor the first three scenarios. If the LN-DLA 
are to be explained through the young age of the absorb- 
ing galaxy, then one would expect the LN-DLA to fill the 
parameter space beneath the plateau instead of compris- 
ing a bimodal distribution. Regarding the second option, 
it would be difficult to explain why DLA with the same 
metallicity exhibit such different N/a values. Finally, the 
third avenue is promising but current models (e.g. Heger 
& Woosley 2002) of Pop III nucleosynthesis appear to un- 
dcrpredict the N/a levels observed in the damped systems. 

Therefore, we favor a scenario where the LN-DLA were 
enriched with an initial starburst with suppressed forma- 
tion of intermediate-mass stars. We found that a starburst 
characterized by a power-law IMF with slope 0.10 or a 



mass cut of w 5 — 8M Q would successfully reproduce the 
observed LN-DLA values. This scenario allows for the bi- 
modal DLA distribution of N/a without requiring that we 
have observed the damped systems at a special moment 
in time. Furthermore, an initial burst of star formation in 
a low metallicity environment may suppress the formation 
of lower mass stars (e.g. Abel, Bryan, & Norman 2000; 
Bromm ct al. 2001) and naturally lend to a bimodal dis- 
tribution. 

Before concluding, we wish to briefly comment on cor- 
relations between N/a and other physical properties of 
the damped Lya systems. In terms of kinematic charac- 
teristics, two of the LN-DLA systems (Q1946+75, z=2.8; 
Q2348-14, z=2.3) have among the simplest kinematics of 
any damped Lya system. This may suggest these sys- 
tems have particularly small mass and/or relatively quies- 
cent star formation. We note, however, that several other 
systems with higher N/a ratios also exhibit very simple 
kinematics. The only other trend we can identify is an 
increase in N/a with increasing z a bs which is contrary to 
expectation. Perhaps this absence of correlation is further 
evidence that the production of N is largely independent 
of environment and star formation history as suggested by 
the uniformity of the emission line sample. 
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